Local and average behavior in inhomogeneous superdiffusive media 
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We consider a random walk on one-dimensional inhomogeneous graphs built from Cantor fractals. 
Our study is motivated by recent experiments that demonstrated superdiffusion of light in complex 
disordered materials, thereby termed Levy glasses. We introduce a geometric parameter a which 
plays a role analogous to the exponent characterizing the step length distribution in random systems. 
We study the large-time behavior of both local and average observables; for the latter case, we 
distinguish two different types of averages, respectively over the set of all initial sites and over the 
scattering sites only. The "single long jump approximation" is applied to analytically determine the 
different asymptotic behaviours as a function of a and to understand their origin. We also discuss 
the possibility that the root of the mean square displacement and the characteristic length of the 
walker distribution may grow according to different power laws; this anomalous behaviour is typical 
of processes characterized by Levy statistics and here, in particular, it is shown to influence average 
quantities. 



PACS numbers: 



I. INTRODUCTION 



The laws of Brownian motion crucially relies on the 
hypothesis that the steps for the diffusing particle are 
small (with finite variance) and uncorrelated. Whenever 
these assumptions are violated, the standard diffusion 

a'cture breaks down and anomalous phenomena emerge 

In particular, transport processes where the step length 
distribution has a diverging variance have been theoret- 
ically studied in detail. Among those, one of the most 
interesting is the so-called Levy walks [3, Q; in which 
particles perform independent steps I at constant veloc- 
ity, with a distribution A(/) following an algebraic tail 
~ ^-(i+a)^ Such a distribution is said to be heavy-tailed 
and it is known to have a diverging variance for a < 2. 
Since transport is thereby dominated by very long steps, 
the mean square displacement increases faster than lin- 
early with time, hence the name superdiffusion. 

Among the many possible experimental applications, 
our work is motivated by the recent realization of ma- 
terials termed Levy glasses, where light rays propagate 
through an assembly of transparent spheres embedded 
in a scattering medium 0, Q- If the diameter of the 
spheres is power-law distributed, light can indeed per- 
form anomalous diffusion. Owing to their ease of fabri- 
cation and tunability, such a novel material offers an un- 
precedented opportunity to study anomalous transport 
processes in a systematic and controllable way. 

An important feature of the experimental samples is 
that the walk is correlated: light that has just crossed a 
large glass microsphere has a higher probability of being 
backscattered at the following step and thus to perform 
a jump of roughly the same length. While the case of 



uncorrelated jumps is well understood [7|, the correla- 
tion effects, that are expected to deeply influence the 
diffusion properties [a], are still to be characterized. To 
this aim, quenched Levy processes have been studied on 
one dimensional systems [9|, |lO| . More recently, different 
aspects regarding the scaling properties of random-walk 
distributions, the relations between the dynamical ex- 
ponents and the different average pro cedures have been 
discussed in a common framework [ll| . 

In order to get a deeper insight on the effect of 
step-length correlations, a class of deterministic, one- 
dimensional models called Cantor graphs has been intro- 
duced [12|. Random walks on these structures perform 
correlated long jumps induced by the underlying frac- 
tal topology. As the latter is generated by deterministic 
rules, diffusion properties can be studied in a simpler way 
than in the random case. Here we extended to this de- 



terministic topologies some of the results proved in 11 1 
for random structures. In particular, we introduce a ge- 
ometric parameter a which plays the same role as the 
exponent characterizing the step length distribution A(/) 
in random systems. Three kinds of statistical averaging 
are introduced: (i) a local one, namely the average of 
all trajectories starting from a given initial site (whose 
asymptotic behaviour being expected to be independent 
of the site choice) ; (ii) an averaging over all possible ini- 
tial sites of the graph and (iii) averages where only scat- 
tering sites are considered as initial conditions. In the 
random case, the differences between such averaging pro- 
cedures have been evidenced in [3, |il[ . On determinisitc 
structures, average procedures (i) and (ii) have been dis- 
cussed in [i^ . Here we complete the picture studying the 
effect of averaging over scattering sites. We evidence that 
the behavior of the mean square displacement is similar 



to random case; while the probability density displays 
a more complex structure being given by a non-trivial 
time-dependent superposition of step functions. We re- 
mark that, in experiments fsl, 'f|, light enters the sample 
with a scattering event and averages of type (iii) are the 
most physically sensible quantities to compare with. In- 
terestingly, in Levy processes the root of the mean square 
displacement and the characteristic length of the distri- 
butions may grow according to different power laws; here, 
in particular, this strongly anomalous behavior charac- 
terizes average quantities. 

The paper is organized as follows. In the next section 
we introduce the directed Cantor graphs and we define 
a simple random walk on these structures. We then dis- 
cuss the relevant physical quantites and the average pro- 
cedures. Section 4 is devoted to discuss the scaling hy- 
pothesis and the single long jump approximation, which 
allows to evaluate the tails of the density distribution. 
Finally, in Section 5 we discuss our results evidencing 
differences among averaging procedures. In particular 
for the case (i) and (ii) we review the results presented 
in |12| presenting new numerical data, while for the case 
of average over scattering sites, which has not been dis- 
cussed so far for the deterministic structures, we provide 
a sketch of the derivation of the dynamical exponents 
within the single long jump approximation and we com- 
pare the result with numerical simulations. In general 
the new simulations evidence that the probability den- 
sity presents different asymptotic behaviors depending 
on the average procedures, clearly supporting both the 
scaling hypothesis and the single long jump approach. 
Moreover the asymptotic behavior of the mean square 
displacement have been tested for a wider range of a's. 



II. RANDOM WALKS ON CANTOR GRAPHS 




FIG. 1: The generation G = 3 of the graph with n„ = 3, 4, 5, 
in the Ur = 2 case. 




FIG. 2: Three examples of Cantor graphs of generation G 
with rir = 2, 3, 4 and n^ ~ 3, 4, 5. 



lowed to move only in one direction (in this case the move 
is performed with probability one, i.e. ballistic motion). 
All the links have unitary length and are crossed at con- 
stant velocity v. The number of bidirectional sites Nb 
and unidirectional sites N^ present at generation G is 
given by: 



Nb = 2nf, 



Nu^ 



(Uu - l)?!^ - [Ur - l)n^ 



-1 (1) 



In paper [12l| we have introduced a class of graphs, de- 
noted as directed Cantor graphs, as a one-dimensional, 
deterministic counterpart of the geometric structure of 
the Levy glass materials mentioned above f5|. Indeed, 
random walks on the Cantor graphs display a superdif- 
fusive Levy like motion analogous to the the one of light 
in such inhomogeneous glassy material. 

The class of graphs we will consider is defined by two 
parameters, denoted as n^ and n„, and describing the 
growth of the fractal from generation G — 1 to generation 
G. In particular the fractal Q of generation G is built 
connecting n^ fractals of generation G — 1 by n^ — 1 uni- 



,G-1 



as shown in 



directional bubbles of length Lq 
Figures [1] and [2] 

A simple random walk 13] is naturally defined on these 
structures: undirected links connect sites in both direc- 
tions while directed links have to be crossed only in the 
prescribed way. A site is called bidirectional if the walker 
placed on that site can move in both directions (and the 
two possible moves are performed with probability 1/2). 
A site is instead called unidirectional if the walker is al- 



where we use the convention of counting only once the 
couple of sites of the bubbles at the same distance from 
the origin. The total number of sites at generation G is 
hence Nq = Nb + N^- In the following, an important 
parameter for the description of the structure will be the 
ratio a — log(n,.)/log(n„). For a < 1, the graphs are 
called slim, as the fraction of bidirectional sites vanishes 
in the thermodinamic limit, i.e. limG->oo Nb/Nc = 0. 
On the contrary, for a > 1 graphs are called fat, since 
a finite fraction of bidirectional sites is present, and 
limc^oo Nb/Na > 0. 



III. PHYSICAL QUANTITIES AND AVERAGES 

Let Pi (r, t) be the probability of arriving at distance r 
starting from i in t steps. In general, Pi{r,t) depends on 
the starting site i. However for large enough space scales 
(i.e. r much larger than the distance between i and an 
origin i = 1) the asymptotic properties of Pi(r,t) are 
expected to be site independent and hence to describe a 



property of the whole graph, i.e. Pi{r,t) ~ Pi{r,t). In 
particular, the asymptotic behaviour of the mean square 
displacement 



with a convergence in probability 



{rUt))= / P,{r,tydr 



(2) 



should be independent of i. The integration cutoff in ([2]) 
is provided by the fact that the walker covers at most a 
distance vt in a time t (v — 1 in the following). Another 
important quantity whose asymptotic behaviour depends 
only on the large scale topology of the structure is the re- 
sistivity Ri (r) i.e the number of bidirectional sites whose 
distance from i is smaller than r. 

On inhomogeneous structures, average and local prop- 
erties are in general different [IJ] , and in structures char- 
acterized by long tails, different averaging procedures can 
yield different results as well p, lll|. In particular, on 
Cantor graphs one can distinguish between the average 
over all starting sites, i.e.: 



P{r, t) ^ lim 

G—^oo 



No 



(3) 



and the averages on processes beginning with a scattering 
event: 



Pir,t) 



lim 






(4) 



where Gb is the set of the bidirectional sites belonging to 
graphs of generation G. The same averaging procedures 
can be introduced also for different quantities. One can 
consider the average resistivity R{r) and R{r) and, as in 
^ , one can also define the average mean square distances 
(f2(i)) and (f2(i)). 

Since the resistivity can be evaluated by simply count- 
ing the number of sites in a given generation, one obtain 
the following asymptotic behaviours: 



i?i(r) r^ R{r 





(5) 



(6) 



These results represent the deterministic analog of the 
expressions for the resistivity obtained in a random sam- 
ple in [10]. 



IV. THE SCALING HYPOTHESIS AND THE 
SINGLE LONG JUMP APPROXIMATION 

The most general scaling hypothesis for the probabili- 
ties Pi{r,t) is: 



P,ir,t) ^ £^\t)hir/h{t)) + g,ir,t) 



(7) 



lim 



\Pi{r,t)-e^\t)h{r/h{t))\dr^Q (8) 



The leading contribution to Pi{r,t) is hence 
^r^(0/i(''/^i(^)) which is significantly different from 
zero only for r < ii{t). The subleading term gi{r,t), 
with limt^oc J \giir,t)\dr = describes the behavior 
at larger distances, i.e. ii{t) -^ r < t. Notice that, if 
5i (r, t) does not vanish rapidly enough, it can never- 
theless provide important contributions to {r1{t)). The 
same scaling ansatz should be valid also for the average 
probabilities P{r, t) and P(r, t) by introducing suitable 
averaged scaling length and scaling functions f{r/£), 
f{r/i), and suitable averaged long distance corrections 

In [l3, ll3 it has been proved that the growth of the 
characteristic length can be directly related to the growth 
of the resistance. We define the exponent describing the 
growth of the correlation length as: 



iW 



dj2 



(9) 



so that in analogy with standard definition of random 
walks [ill, we get Pi{0,t) ~ t''^^/^. Then using the 
scaling relations proved in [1^, Hal we obtain 



Mr) 



^2/d,,-l 



(10) 



Analogous relations hold for the average quantities, in 
terms of the average scaling lengths. Introducing the 
known results for the resistivity (|5l6p . one obtains the 
following behaviors for the scaling lenghts: 



ii+° , if a < 1 



ei{t)^m 



m 



a 



if a> 1 



t if a < 1 



t2 if a > 1 



(11) 



(12) 



Let us now discuss the behaviour of the mean square 
displacements. When only lengths of order r < i{t) pro- 
vide significant contributions to the integral ([2]), the stan- 
dard relation (r^(i)) ^ £'^{t) holds and the asymptotic be- 
haviour coincides with those given by (|11I12[) . However 
it is known that, in presence of long tailed distributions, 
anomalies with respect to this behaviour can be present. 
In particular it has been evidenced in [^, [ll| that a key 
role is played by long jumps, leading the walker to a dis- 
tance r ^ £{t). In the random case, consideration of 
a single loiig jump actually accounts for the asymptotic 
behaviour [li|. Indeed, these processes can give rise to 
two different types of corrections to P{r,t). First, they 
can produce a zero- measure function g{r,t), providing 
a significant contribution to (r^(i)); second, the scaling 
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FIG. 3; Dynamical scaling of the probability for initial site 
i = 1 on the graph obtained with rir = 2 and n^ = 3. Note 
the fast decay of the scaling function. The stepwise structure 
of the scaling function is due to the fractality of the graph. 



FIG. 4: Growth of the mean square displacements for ini- 
tial site i = 1. Power-law fits give exponents in very good 
agreement with the theoretical values, in equation p3p . 



function /(x) can feature a long tail, breaking the pro- 
portionality between (r^(i)) and i'^{t). Here we will ev- 
idence that these anomalies, originating from the single 
long jump, are also present in the deterministic graphs 
and they are deeply influenced by the averaging proce- 
dures. 



in Figure HI the continuous lines represent the expected 
behaviors (IT^ . 



B. Averages over all sites 



V. RESULTS 



A. Local behaviour 



Let us first consider the local properties. In this situ- 
ation we can focus on processes starting from the origin 
j = 1 of the graph. Indeed, if r is much larger than 
the distance between i and 1, we expect that Pi{r,t) 
should behave as Pi (r, t) and hence the asymptotic be- 
haviour should be the same for any starting point. When 
starting from the origin, in a time t the walker typically 
covers a distance ii{t) and, in the deterministic graphs, 
within such distance there are no bubbles of length larger 
than ii(t). Therefore, long jumps do not occur and 
(r?(i)) ~ ejit). Therefore [l|| 



(13) 



Figure [3] reports, an example the probability density 
Pi (r, t) obtained by a Montecarlo simulation. The data 
provide a clear evidence that the scaling function presents 
a fast decay confirming our hypothesis of no long jump. 
As explained in [12] the fractal structure give rise to log- 
periodic oscillations which can be discarded considering 
peculiar sequence of times (in this case t = 6*^) , However, 
such oscillations do not change the general framework of 
the scaling hypothesis. The growth of (rj(i)) is plotted 




When averaging over all the sites, it has been evidenced 
[13 | that for a < 1 the motion is always ballistic, while 
for a > 1 the situation is much more complex, since the 
walker can perfom a single jump much larger than £(t). 
Tipically, such a long jump occurs at the first step be- 
cause, with a random choice of the starting point, the 
probability of belonging to a large bubble is much larger 
at i = than during the rest of the evolution. In this sit- 
uation, one can estimate the zero measure correction to 
P{r,t) obtaining g{r,t) ~ i~"+^(5(r — i), i.e. a peak asso- 
ciated with the ballistic motion of the particle, weighted 
by a factor t~°'~^^ representing the probability of belong- 
ing to a bubble larger than t at the initial time. The be- 
havior of P{r, t) is illustrated in the simulations of figure 
[5] evidencing the presence of a scaling regime for r < £{t) 
and of balistic peaks at large r whose height evolves as 
^-Q+i ^(iased line). Even if ^(r, i) provides a subleading 
contribution to P{r, t), the integral in ^ is dominated by 
g{r,t) for 1 < a < 2 and by f{r/£{t)) only when a > 2. 
Therefore the behaviour of the mean square displacement 



(r'it)) 



i^-", if 1 < a < 2 



if a>2 



(14) 



while the motion is purely ballistic for a < 1. Figure [5] 
shows that the predicted exponents P^ are well verified. 





FIG. 5: Dynamical scaling of P{r, t) in the average case. The 
data refer to the case n^ = 7 n„ = 4. Ballistic peaks at r = f 
scale as f^~" (dashed line), while for r < £{t) the scaling 
hypothesis is well satisfied. 



FIG. 7: The dynamical scaling of P{r, t) in the case of aver- 
ages over bidirectional sites. Here we plot the case n^ = 4 
and n„ = 3. For r < l{t) the scaling hypothesis (0 is very 
well verified with t.(t) growing as indicate in (|lip . For larger 
r the behaviour is much more complicated since the tails are 
composed by a superpositions of step functions. 
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FIG. 6: Growth of the mean square displacements when aver- 
age over scattering over the whole sample is considered. We 
compare the numerical results with the theoretical preditions 
(fT4|l in both the regimes 1 < a < 2 and q > 2. 



C. Averages over bidirectional sites 

In [ll| it has been shown that, for random systems, 
averages over bidirectional scattering sites can provide 
different results with respect to averages over all starting 
points and a scaling approach has been discussed, based 
on the single long jump approximation. Here we intro- 
duce an analogous argument for the deterministic graph. 

In the case of average over bidirectional sites, the single 
long jump does not occur necessarily at the first step. In 
particular for the deterministic graph the probability of 
performing a jump of length Lk = n^^ £{t) in a time t is 
N{t)n~'' where N{t) is the number of bidirectional sites 
visited by the walker in a time t and n~ is the probability 
that a bidirectional site belongs to a bubble of length Lk- 
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FIG. 8: Growth of the mean square displacement when av- 
eraged over bidirectional sites is consider. The results are 
compared with the theoretical predictions (|18|) given by the 
continuous lines. Different values of n^ and n„ are considered, 
showing that the predictions (|18|) are verified in the different 
regimes of q's. 



Discarding this long jump, the distance crossed by the 
walker in a time t is of order £{t), and therefore, according 
to the behaviour of the resistivity described by equations 
B N{t) - £(i)" for a < 1 and N{t) ~ e(t) for a > 1. 
The main difference with respect to the random case is 
that now the only possible lengths of the long jumps are 
given by the sizes Lk of the bubbles in the fractal. Hence, 
for i{t) < r < t, P{r, t) is a complex step function where 
both integers n„ and rir plays a non trivial role. However, 
the contribution to the mean square displacement can be 



evaluated as follows: 



m\ E ■ 

ye(t)<Lk<t 



^2k 



+ r 



Lk>t 



(15) 



where the first sum is related to the bubbles of length 
Lk i^t) < Lk < t) providing a contribution to the mean 
square displacements of order L\ = n^ , while the sec- 
ond sum comes from the bubbles of length larger than i, 
providing a contribution t^ . Expanding equation [TS] for 
large times one obtains the asymptotic behaviours 



2fe 



m\ E ^+^^Ei 

\i(t)<Lk<t '' 
for a < 1 and 



2 + 2e, + Q 
t 1 + ° 



(16) 



Lk>t 



Equations (|18p extend the results of the random case 
to the deterministic topology described by the directed 
Cantor graphs. Clearly the complex shape of the scal- 
ing function (|15p determines the presence of logperiodic 
oscillations superimposed to (1181) . which is a a typical 
behaviour of fractal structures [131 • Figure (O evidences 
that the dymanical scaling is well verified for r < ^(i), 
while at larger distances P(r, t) is characterized by a su- 
perposition of slowly decaying step functions as predicted 
by (ITS]) . Figure ^ evidences by means of Montecarlo 
simulations that Equations (TTSl) are well verified for large 
times in the whole range of a's. 



N{t) 



E 



-,2k 



\t{t)<Lk<t 



Lk>t ^ 



5/2-Q 



(17) 



for a > 1. The first expression is always dominant with 
respect (.{t), while the second expression becomes sub- 
leading for a > 3/2. The overall behaviour of the mean 
square displacement is summarized as follows: 



if < a < 1 

if 1 < a < 3/2 (18) 

if 3/2 < a 
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